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Abstract.  This  paper  is  concerned  with  the  mathematical  analysis  of  the  time-domain  Maxwell 
equations  in  a  three-dimensional  cavity.  An  exact  transparent  boundary  condition  (TBC)  is  devel¬ 
oped  to  reformulate  the  open  cavity  scattering  problem  in  an  unbounded  domain  equivalently  into  an 
initial-boundary  value  problem  in  a  bounded  domain.  The  well-posedness  and  stability  are  studied 
for  the  reduced  problem.  Moreover,  an  a  priori  estimate  is  established  for  the  electric  field  with  a 
minimum  regularity  requirement  for  the  data. 
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1.  Introduction.  This  paper  is  concerned  with  the  mathematical  analysis  of 
an  electromagnetic  cavity  scattering  problem  governed  by  the  time-domain  Maxwell 
equations.  An  open  cavity  can  be  broadly  defined  as  a  local  perturbation  of  an  infinite 
plane  surface.  The  cavity  scattering  problems  have  significant  industry  and  military 
applications  such  as  the  design  of  cavity-backed  conformal  antennas  and  the  deliberate 
control  in  the  form  of  enhancement  or  reduction  of  radar  cross  section. 

The  time-harmonic  problems  have  been  widely  investigated  by  numerous  re¬ 
searchers  in  the  engineering  and  mathematical  communities  [2-4,7,22,24,30].  A 
large  amount  of  information  is  available  regarding  their  solutions  for  both  the  two- 
dimensional  Helmholtz  and  the  three-dimensional  Maxwell  equations  including  the 
overfilled  cavity  problems  [7,8,27,29,41].  We  refer  to  [23]  for  a  good  introduction 
to  the  cavity  scattering  problem.  One  may  consult  [14,15,31,32]  for  recent  accounts 
of  finite  element  methods  and  integral  equation  methods  for  general  electromagnetic 
scattering  problems. 

The  time-domain  problems  have  attracted  much  attention  due  to  their  capability 
of  capturing  wide-band  signals  and  modeling  more  general  material  and  nonlinearity 
[9,21,26,39,40,42].  However,  rigorous  mathematical  analysis  is  very  rare,  especially 
for  time-domain  three-dimensional  Maxwell  equations.  The  transient  cavity  scattering 
problems  have  been  mathematically  studied  in  [19,20,34,36,37],  where  the  focus  was 
largely  on  (i)  temporal  discretization  by  e.g.,  the  Newmark  scheme;  (ii)  the  reduction 
of  the  resulted  system  via  frequency-domain  TBC,  and  (iii)  the  analysis  of  the  finite- 
element  method  for  the  reduced  problem.  To  the  best  of  our  knowledge,  the  theoretical 
analysis  of  the  time-domain  Maxwell  cavity  problem  itself  was  left  undone  and  is  still 
lacking. 

In  this  work,  we  aim  to  analyze  the  transient  electromagnetic  wave  scattering  by 
an  open  cavity  which  is  embedded  in  a  perfectly  conducting  infinite  ground  plane. 
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The  cavity  is  filled  with  some  inhomogeneous  medium  which  may  even  protrude  from 
its  opening  to  above  the  ground  plane,  while  the  upper  half  space  outside  of  the  cavity 
is  filled  with  a  homogeneous  medium.  This  problem  is  quite  challenging  due  to  the 
unbounded  nature  of  the  domain.  In  the  frequency  domain,  various  approaches  have 
been  investigated  such  as  developing  TBCs  and  designing  perfectly  matched  layer 
techniques  [1,5,10,12,17,18,25,28,29].  In  contrast,  much  less  study  is  devoted  to 
address  this  issue  in  the  time  domain.  In  [11],  the  well-posedness  and  stability  were 
shown  for  the  electromagnetic  obstacle  scattering  problem  by  using  a  time-domain 
TBC  on  the  sphere.  Here,  we  present  a  new  time-domain  transparent  boundary 
condition  on  the  hemisphere  emcompassing  the  cavity  and  prove  well-posedness  and 
stability  of  the  underlying  scattering  problem.  The  proofs  are  based  on  examining  the 
well-posedness  of  the  time-harmonic  Maxwell  equations  with  complex  wavenumbers 
and  applying  the  abstract  inversion  theorem  of  the  Laplace  transform.  Moreover,  an  a 
priori  estimate  featured  with  explicit  dependence  on  time  and  a  minimum  regularity 
requirement  of  the  data,  is  established  for  the  electric  field  by  studying  directly  the 
time-domain  Maxwell  equations. 

The  outline  of  this  paper  is  as  follows.  In  section  2,  we  introduce  the  model 
problem  and  exploit  the  time-domain  TBC  to  reduce  the  underlying  problem  to  an 
initial-boundary  value  problem  in  a  bounded  domain.  In  section  3,  we  analyze  two 
auxiliary  problems  pertinent  to  the  reduced  time-domain  Maxwell  equations  to  pave 
the  way  for  the  analysis  of  the  main  results  in  section  4.  We  study  in  section  4  the 
well-posedness  and  stability  of  the  reduced  Maxwell  equations,  and  derive  an  a  prior 
estimate  with  a  minimum  requirement  for  the  regularity  of  the  data.  The  paper  is 
concluded  with  some  remarks  and  directions  for  future  work  in  section  5.  To  avoid 
distraction  from  the  main  results,  we  collect  in  the  appendices  some  necessary  nota¬ 
tion  and  useful  results  on  the  Laplace  transform,  spherical  harmonics,  and  functional 
spaces. 

2.  Formulation  and  reduction  of  the  problem.  In  this  section,  we  intro¬ 
duce  the  mathematical  model  of  interest,  and  exploit  the  exact  TBC  to  reduce  the 
unbounded  domain  to  a  bounded  one. 

2.1.  A  model  problem.  We  describe  the  setting  of  a  cavity  problem  and  de¬ 
fine  some  necessary  notations.  Denote  by  D  the  cavity  embedded  in  the  perfectly 
electrically  conducting  infinite  ground  plane  Tg.  Let  S  =  dDf lK_,  the  part  of  cavity 
walls  below  the  ground,  be  Lipschitz  continuous  and  perfectly  electrically  conducting. 
The  cavity  is  filled  with  some  inhomogeneous  medium,  which  may  protrude  from  its 
opening  to  above  the  ground  plane.  The  medium  is  characterized  by  the  dielectric 
permittivty  £  and  the  magnetic  permeability  p,  which  satisfy 

0  A  £min  —  ^max  A  OO,  0  <C  p min  —  P  —  Mmax  A  OO. 

Here  £min,  fmax,  Mmim  and  p-max  are  constants.  Therefore,  the  analysis  in  this  work 
includes  both  planar  and  overfilled  cavities.  Let  and  rjj  be  the  half-ball  and 
hemisphere  above  the  ground  plane,  where  the  radius  R  is  large  enough  to  completely 
contain  the  possibly  overfilled  portion  of  the  cavity.  The  unbounded  region  \  B ^  is 
filled  with  some  homogeneous  material  with  a  constant  permittivty  £o  and  a  constant 
permeability  po-  Throughout  this  paper,  we  assume  for  simplicity  in  exposition  that 
£o  =  Po  =  1-  Finally,  we  denote  by  fl  =  B ^  U  D  the  bounded  domain  in  which  our 
reduced  initial-boundary  value  problem  will  be  formulated. 

It  is  evident  that  the  problem  geometry  is  not  only  applicable  to  the  cavity  prob¬ 
lem  but  also  to  a  broader  class  of  scattering  problems  where  the  interface  or  the 
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Fig.  2.1.  ^4  schematic  diagram  of  the  cavity  problem  geometry. 


boundary  is  modeled  as  a  local  perturbation  of  an  infinite  plane. 

Consider  the  system  of  time-domain  Maxwell  equations  in  U  D  for  t  >  0: 


V  x  E(r,  t)  +  iidtH (r,  t)  =  0, 

V  x  H(r,t)  —  edtE(r,t )  =  J(r,t), 


where  r  =  (. x,y,z )  G  M3,  E  is  the  electric  field,  H  is  the  magnetic  field,  and  J  is 
the  electric  current  density  which  is  assumed  to  be  compactly  supported  in  D.  The 
system  is  constrained  by  the  initial  conditions: 

(2.2)  E\t—o  =  Eo,  H  |t=0  =  H0  in  R^UD, 

where  E o  and  H 0  are  also  assumed  to  be  compactly  supported  in  D.  We  consider 
the  perfectly  electrical  conducting  boundary  condition  on  the  ground  plane  and  cavity 
wall: 


(2.3)  n  x  E  =  0  on  Ts  U  5,  t  >  0, 

where  n  is  the  unit  outward  normal  vector  on  Tff  U  S.  In  addition,  we  impose  the 
Silver-Miiller  radiation  condition: 

(2.4)  f  x  (dtE  x  r)  +  r  x  dtH  =  o(|r|_1),  as  |r|  — >  oo,  t  >  0, 
where  f  =  r/\r\. 

The  purpose  of  this  paper  is  to  study  the  well-posedness  and  establish  the  stability 
for  the  time-domain  electromagnetic  cavity  scattering  problem  (2.1)-(2.4).  Hereafter, 
the  expression  “a  <  b”  stands  for  a  <  Cb ,  where  C  is  a  generic  positive  constant 
independent  of  any  function  and  important  parameters,  which  are  clear  from  the 
context. 

2.2.  Transparent  boundary  condition.  We  introduce  a  TBC  to  reformulate 
the  electromagnetic  wave  propagation  problem  from  an  open  cavity  into  an  equivalent 
initial-boundary  value  problem  in  a  bounded  domain.  The  essential  idea  is  to  design 
a  boundary  operator  which  maps  the  tangential  component  of  the  electric  field  to  the 
tangential  trace  of  the  magnetic  field. 
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More  precisely,  we  consider  the  reduced  initial-boundary  value  problem: 


Vx£  +  pd tH  =  0,  V  x  H  -  edtE  =  J 

in  Q,  t  >  0, 

E  |t=o  =  Eq,  H  |t=o  =  i?o 

in  fl, 

n  x  E  =  0 

on  S,  t  >  0, 

8T[Et8\  =H  xn 

on  T+,  t  >  0 

where  Er+  is  the  tangential  trace  of  E  on  T^,  and  88  is  the  time-domain  electric-to- 
magentic  (EtM)  Calderon  operator,  as  the  counterpart  of  the  time-harmonic  setting 
with  e.g.,  the  spherical  boundary  (cf.  [15,31]). 

In  what  follows,  we  derive  the  formulation  of  the  operator  88  and  analyze  its 
important  properties.  Equivalently,  we  aim  to  prove  the  well-posedness  and  stability 
of  the  reduced  problem  (2.5).  In  particular,  a  priori  estimates  are  established  with 
minimum  requirements  for  the  regularity  of  the  data.  The  proofs  are  based  on  the 
abstract  inversion  theorem  of  the  Laplace  transform  and  the  a  priori  estimates  for 
the  time-harmonic  Maxwell  equations  with  a  complex  wavenumber.  These  will  be  the 
main  topics  of  the  forthcoming  sections. 

Since  J  is  supported  in  D  and  £o  =  Mo  =  1,  the  Maxwell  equations  (2.1)  become 

(2.6)  \7xE  +  dtH  =  0,  V  x  H  -  dtE  =  0  in  fle  :=  M3.  \  B+,  t  >  0. 

Let  E(r,  s )  =  88(E)  and  H(r,  s )  =  88(H)  be  the  Laplace  transforms  of  E(r,  t)  and 
H(r,t)  with  respect  to  t,  respectively,  where  the  complex  variable 

(2.7)  s  =  Si  +  is2  with  si,S2  €  K,  si  >  0,  i  =  y/—  1 , 
as  seen  in  Appendix  A.  Recall  that 

88(dtE)  =  sE  -  E0,  88(dtH)  =  sH  -  H0. 


Taking  the  Laplace  transform  of  (2.6),  and  noting  that  Eq,Hq  are  supported  in  D, 
we  obtain  the  time-harmonic  Maxwell  equations  with  complex  parameters: 

(2.8)  V  x  E  +  sH  =  0,  V  x  H  —  sE  =  0  in  Oe,  sq  >  0. 


Let  hn\z)  be  the  spherical  Hankel  function  of  the  first  kind  of  order  n  (cf.  [38]). 
We  introduce  the  vector  wave  functions 


MZ(p,0,<p)  =  V  x  (rh£\isp)X™(0,<p)), 
N™(p,0,<p)  =  -s-1VxM™(e,<p), 


where  (p,  9,  tp)  is  the  spherical  coordinates,  and  X™  is  the  rescaled  spherical  harmonic 
function  defined  in  (B.2).  We  refer  to  Appendix  B  for  the  properties  of  spherical 
harmonics  and  the  relevant  calculus  to  be  used  throughout  this  paper. 

The  vector  wave  functions  in  (2.9)  are  the  radiation  solutions  of  (2.8)  in  M3  \  {0} 
(cf.  [31,  Thm.  9.16]): 

(2.10)  V  x  M™(p,9,cp)  +  sN™(p,0,<p)  =  0f  VxN™(p,0,cp)-sM™(p,0,lp)  =  O. 
It  can  be  verified  from  (2.9)  that  the  vector  wave  functions  satisfy 

(2.11)  M™(p,e,<p)  =  h«(iSp)VrX-(0,^)  x  ep 


Transient  Electromagnetic  Scattering 


5 


and 

N™(p,e,ip)  =  -y/n^  +  1)  (hW(isp)+isp(hPnispj)  x™{e,p) 

(2.12) 

Once  again,  we  refer  to  Appendix  B  for  the  notation  and  definition.  A  simple  calcu¬ 
lation  yields 

{eg  x  M™(p,  7t/2,  ip)  =  0  for  \m\  <  n,  m  +  n  =  even,  n  £  N, 

eg  x  N™(p,  7t/2,  (p)  =  0  for  \m\  <  n,  m  +  n  =  odd,  n  £  N. 

Therefore,  in  the  domain  Qe ,  the  solution  of  the  electric  field  E(p,0,(p,s),  which 
satisfies  the  perfectly  electric  conducting  condition  n  x  E  =  0  on  f9,  i.e. ,  eg  x 
E(p,  7t/2,  p,  s)  =  0,  can  be  written  in  the  series: 

odd  even 

(2.13)  E(p,0,ip,s)  =  ^  o%{s)NZ(p,0,<p)+  E 

|m|  <n  |m|<n 

which  is  uniformly  convergent  on  compact  subsets  in  Qe.  The  corresponding  magnetic 
field  H  is  given  by 

odd  even 

(2.14)  H  =  — s_1V  x  E  =  —  E  a™(s)M™(p,0,<p)  +  E  /TO  N^(p,0,<p)- 

|m|<n  |m|<n 

Note  that  in  the  above,  we  compressed  the  summation  notation  as  in  Appendix  B. 

To  deduce  the  explicit  representation  of  the  EtM  Calderon  operator,  we  need  to 
express  Er+  =  —ep  x  (ep  x  E)  and  H  x  ep  on  T^  in  terms  of  the  coefficients  a™  and 
/3™.  From  the  definition  (2.11),  one  verifies  that 

-eP  x  (ep  x  M%(p,0,<p))  =  -^n{n  +  l)h^)(\sp)Y^{01p), 

~eP  x  (ep  x  N™(p,  0,  ip))  =  -  +  ^  (h£\isp)  +  isph^' (isp)j  X™{0,<p), 

and 

ep  x  M™(p,0,ip)  =  yj n(n  +  1)  h ^  (i sp)X™(0,  ip), 

ep  x  N™{p,0,p)  =  -  1'>-  (h£\isp)+isph£)\isp)S)  Y^(0,ip). 

Therefore,  by  (2.13),  the  tangential  component  of  the  electric  field  along  T^  is 
£r+=-  E  ^  (h£\isR)+isRh£y(isR))a™(s)X™(0,ip) 

|m|  <n 

even 

|m|<n 
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and  similarly,  by  (2.14),  the  tangential  trace  of  the  magnetic  field  along  T^  is 


odd 


Hxep=  Y  Vn(n  +  1)  h^(iaR)a^(a)  X™(9,  <p) 


|m|  <n 
even 


+  E  VnisR  ^  (h^(isR)  +  isRh^\isR))  f3™(s)Y™(e,<p). 

|m|  <n 

Consequently,  we  have  the  following  explicit  representation  of  the  frequency  do¬ 
main  EtM  Calderon  operator  SS\  given  any  tangential  component  of  the  electric  field 
along  with  the  expansion: 


(2.15) 


odd  even 

U=  E  <xn  +  E  WY 

|m|<n  |m|<n 

the  tangential  trace  of  the  magnetic  field  along  T)|)  is 

odd  7-,  even 

sn 


(2.i6)  a\u]  =  -  y  — E  (1+r"l(isi?))^^, 

\m\<n  1  +  r«  (1SR) 


|m|  <n 


sR 


where 

(2.17) 


ID 

"  (  }  h£\z)  ' 


We  now  analyze  some  properties  of  the  EtM  Calderon  operator,  and  refer  to 
Appendix  D  for  the  definitions  of  the  function  spaces  to  be  used  hereafter. 

Lemma  2.1.  The  Calderon  operator  £$  :  iT_1/2(curl,  T^)  —>  iL-1/2) div,  T^)  is 
continuous. 

Proof.  For  any  u,v  £  ff~1^2( curl,  T^),  we  can  expand 

odd  even  odd  even 

U=  E  <nXn+  E  V  =  £  <  X™  +  £  V?n  Y™ . 

|m|<n  |m|<n  |m|<n  |m|<n 

Then  we  have  from  (2.16)  that 

odd 


#[«]  =  -  E 


sRu 


In 


\m\<Ln  1  +  )  (1S^)  |m[< 

To  prove  the  lemma,  it  is  required  to  estimate 

odd 


-x™-  E 


{l  +  r{n]  {\sR))u^n 


yrm 


(2.18)  {@[u],v)  =  -  Y 


sR  um  Vm  -  V  (1+r”  m 

^  1  ,  (l)r  u\  ln  1  n  2-^  SR  U2n^2n- 

|  m\<n1+rn>(lsR)  |m|<n 


It  follows  from  the  Cauchy- Schwarz  inequality  that 


\{3§[u),v)\  < 


ST'  y/l  +  n(n  +  1)  I  n|2  I  m  |2  _j_  Y"  I1  +  E\isA)|2  I'^l2 

^  |1  +  r«  W  1  1  1  lJ  +  >A  +  »(»  +  D  W 


1/2 


odd  1  even 

E  77=  7  ,-.>ml2+  E  >/l  +  n(n+l)|« 

,  ...  \/l  +  n(n+l  ,  ,, 

|m| <n  v  v  '  lm.l<n. 


1  1/2 


m  |2 
2n  I 


|m|<n 
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By  Lemma  C.3,  we  have 


and 


|1  +  r 


n\n 

(i) 


!)  i  Di2i,m|2  1  1  +  n(n  +  1)  |  |2Lm, 

7 - 777  \sR\  \uln |  =  ,  .  ==: - m7 - —\sR\  \uin\ 

(isf?)|2  \/l  +  n(n  +  1)  |1  +  r„  (isf?)|2 

1 


< 


\/l  +  n(n  +  1) 


Knl2, 


ll  +  r^^isfl)!2  fe|2  _  /  ~7~~  1  Jl  +  r-^Qs-R)!2  |ug,|2 
^/l  +  n(n  +  l)  1^1 2  V  1  J  1  +  n{n  + 1)  |sf?|2 

<  a/1  +n{n+  1)  \u™n\'2. 


Combining  the  above  estimates  and  using  the  expressions  of  the  norms  in  Appendix 
D  yields 


\{@[u),v)\  ~  IMI  j/— 1/2(curi,r+)  llt;llff-1/2(CUri, r+)> 
which  completes  the  proof.  □ 

Another  important  property  of  the  EtM  Calderon  operator  is  stated  as  follows. 
Lemma  2.2.  It  holds 

(2.19)  Re  (St[u\,u)  >  0,  Vm  G  ff~1/2(curl,  T+). 


Proof.  From  (2.18),  we  obtain 


odd 


(2.20)  -<#[«],«>=  £ 

|m|<n 


si? 


1  +  (isi?) 


m  |2 
Anl 


1  +  rn^isi?)  |2 


|m|<r: 


sR 


x2nl 


By  Lemma  C.l  and  Lemma  C.2, 

SiRe(l  +  r^^isi?))  <  0,  s2Im(l  +  r^(isi?))  <  0. 
Taking  the  real  part  of  (2.20)  gives 


|m|<n 


|1  +  ri1')(isi?)|2 


n  I 


^  (siRe(l  +  ri1)(isi?))  +  s2Im(l  +  ri1)(isi?)))|  m  |2  ^  n 

/  ,  U2D  lU2nl  A  U, 


|m|<r; 


s2R 


which  completes  the  proof.  □ 

With  the  aid  of  the  frequency  domain  EtM  Calderon  operator,  we  obtain  the 
following  TBC  imposed  upon  the  hemisphere  T^  in  the  s-domain: 

(2.21)  £8[Er+\  =  H  x  ep, 

which  maps  the  tangential  component  of  the  electric  field  to  the  tangential  trace  of 
the  magnetic  field.  Taking  the  inverse  Laplace  transform  of  (2.21)  yields  the  TBC  in 
the  time-domain: 


(2.22)  S'[Er+\  =  H  x  ep,  where  S'  :=  Jz?  1  o  St  o  ££ . 


P.  Li,  L.  Wang,  and  A.  Wood 


Equivalently,  we  may  eliminate  the  magnetic  field  and  obtain  an  alternative  transpar¬ 
ent  boundary  condition  in  the  s-domain: 

(2.23)  s-1(V  xE)xn  +  @[Er+]  =  0  on  T+. 

Correspondingly,  by  taking  the  inverse  Laplace  transform  of  (2.23),  we  may  derive  an 
alternative  TBC  in  the  time  domain: 

(2.24)  (V  x  E)  x  n  +  e^’[Er+]  =  0  on  where  o  sC%  o  Jz?. 

3.  Analysis  of  two  auxiliary  problems.  In  this  section,  we  make  necessary 
preparations  for  the  proof  of  the  main  results  by  considering  two  auxiliary  problems 
pertinent  to  (2.5). 

3.1.  Time-harmonic  Maxwell  equations  with  a  complex  wavenumber. 

This  subsection  is  devoted  to  the  mathematical  study  of  a  time-harmonic  Maxwell 
scattering  problem  with  a  complex  wavenumber,  which  may  be  viewed  as  a  frequency 
version  of  the  initial-boundary  value  problem  of  the  Maxwell  equations  under  the 
Laplace  transform. 

Consider  the  auxiliary  boundary  value  problem: 

{V  x  ((s/z)-1V  x  u)  +  ssu  =  j  in  Cl, 
n  x  u  =  0  on  S, 

s_1(V  x  u)  x  n  +  J[«r]  =  0  on 

where  s  =  si  +  is2  with  si,S2  £  K,  si  >  0,  and  the  applied  current  density  j  is 
assumed  to  be  supported  in  D  C  Cl. 

By  multiplying  a  test  function  v  £  Li's  (curl,  Cl),  which  is  defined  in  Appendix 
D,  and  integrating  by  parts,  we  arrive  at  the  variational  formulation  of  (3.1):  find 
u  £  if  5  (curl,  Cl)  such  that 

(3.2)  aTH(u,u)  =  /  j  vdr  for  all  v  £  if 5 (curl,  Cl), 

J  n 

where  the  sesquilinear  form 

(3.3)  a-r}i(u,v)  =  /  (s/i)_1(Vxu)'(Vxi)dr+  /  seu  ■  v  dr  +  {8§[uT+],vT+) . 

Jn  J  n  r  r 

Theorem  3.1.  The  variational  problem  (3.2)  has  a  unique  solution  u  £  Hs(cwrt,Cl), 
such  that 


II V  X  m||i,2(q)  +  ||sM||£2(n)  <  s1 1||sj||I/2(0), 

where  s  =  s  1  +  is2  with  si,  S2  £  R.  and  s  1  >  0. 

Proof.  It  suffices  to  show  the  coercivity  of  axH)  since  its  continuity  follows  di¬ 
rectly  from  the  Cauchy-Schwarz  inequality,  Lemma  2.1,  and  Lemma  D.l.  A  simple 
calculation  yields 


=  /  (sp.)  1|Vxu|2dr+  /  se|w|2dr  + 


u 


]>  Mr+)- 


(3.4) 


Oth(m,  u) 
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Taking  the  real  part  of  (3.4)  and  using  Lemma  2.2,  we  get 


(3.5) 


Re{aTH(u,u)}  >  r^(||V  x  u\\l2{q) 


sm 


L2(Q) 


)• 


It  follows  from  the  Lax-Milgram  lemma  that  the  variational  problem  (3.2)  has  a  unique 
solution  u  £  Hg (curl,  O).  Moreover,  we  have  from  (3.2)  that 


(3.6)  |uth(m,m)|  <  |s|  1\ 

Combining  (3.5)-(3.6)  leads  to 


i2(o)lls'ullz,2(n)- 


IIV  x  «|||,3(n) 


ISUllracoi  5,  1 


i2(o) 


1  llsJlli2(Q)llSMlli2(0)) 


which  completes  the  proof  after  applying  the  Caucliy-Schwarz  inequality.  □ 

It  is  noteworthy  that  Theorem  3.1  gives  the  stability  estimate  with  an  explicit 
dependence  on  the  complex  wavenumber  s  correponding  to  a  lossy  medium  with 
Si  >  0.  However,  it  is  challenging  to  obtain  a  similar  estimate  for  a  lossless  medium 
with  si  =  0. 

3.2.  A  second  auxiliary  problem.  Before  considering  the  reduced  problem 
(2.5),  we  need  to  study  an  auxiliary  initial-boundary  value  problem  and  establish  its 
wcll-posedness  and  stability.  Consider  the  system  of  time-domain  Maxwell  equations: 


V  x  U  +  ndtV  =  0,  V  x  V  -  sdtU  =  0 
(3.7)  |  n  x  U  =  0 

U\t=o  =  E0,  V|t=o  =  H0 


in  f 2,  t  >  0, 
on  SUT)j,  t  >  0, 
in  fl, 


where  Eq.  Hq  are  assumed  to  be  compactly  supported  in  D  as  before. 

Let  U  =  jSf(f/)  and  V  =  Jf(V).  Taking  the  Laplace  transform  of  (3.7),  we 
obtain  the  boundary  value  problem: 


(3.8) 


,Vx((s/j)-1Vxl/)+sE[(=j  in  fi 
n  x  U  =  0 


on  s u rj, 


where  the  current  density  function 

j  =  eE0  +  s^V  x  H0. 

The  variational  formulation  for  (3.8)  is  to  find  U  £  ifo(curl,  f2)  such  that 

(3.9)  a,Ap(U,v)  =  /  j  ■  vdr  for  all  v  £  if o (curl, fl), 

Jn 

where  the  sesquilinear  form 

(3.10)  oap(m,  v)  =  /  (s/j)^1(Vx«)'(Vxn)dr+  /  seu-v&r. 

Jci  Jn 

Following  the  same  proof  as  in  Theorem  3.1,  we  can  show  the  well-posedness  of 
the  variational  problem  (3.9)  and  its  stability,  as  stated  below. 
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Lemma  3.2.  The  variational  problem  (3.9)  has  a  unique  solution  U  €  i?o(curl,  11). 
Moreover,  we  have  the  estimate 

II V  x  +  HsfrH^n)  <  1|s|||^o||x,2(n)  +  «i  1||V  x  i?o||i2(n)- 

Theorem  3.3.  The  auxiliary  problem  (3.7)  has  a  unique  solution  (U,V),  which 
satisfies  the  stability  estimates 

ll^rlli2(n)  +  ll''/'lli2(n)  <  ll-Eolk2(n)  +  II-Ho||z,2(Q)> 

\\dtU\\L2[n)  +  H^VIU^nj  <  ||V  x  E0\\LHn)  +  ||V  x  H0|li2(n), 

\\d?U\\LHn)  +  \\d2tV\\L*{il)  <  ||V  x  (V  x  So)IL»(n)  +  ||V  x  (V  x 


Proof.  Let  U  =  J£(U)  and  V  =  Jf(V)  as  before.  Taking  the  Laplace  transform 
of  (3.7)  leads  to 


(3.11) 


V  x  U  +  spV  =  pH o,  V  x  V  —  seU  =  —e E0 

n  x  U  =  0 


in  12, 

on  snr+. 


It  follows  from  Lemma  3.2  that 

II V  x  U ||i,2(fj)  +  ||sC7||i,2(n)  <  1|s|||£lo||£2(n)  +  «i  1||V  x  iTo||z,2(n)- 

By  (3.11),  we  have 

l|Vx V||i2(n)  +  ||sV'||i2(n)  <  (l+sr'laDllBolU^nj  +  IIHoll^^+sr'llVxHolU2^). 

It  follows  from  Lemma  A. 2  that  U  and  V  are  holomorphic  functions  of  s,  and  the 
inverse  Laplace  transform  of  U  and  V  exist  and  are  supported  in  [0,  oo). 

We  next  prove  the  stability.  Define  the  energy  function 

ei  (t)  =  \\£1/2U(-,t)\\2L2{n)  +  \\v1/2V{-,t)\\2L2{ny 
It  follows  from  (3.7)  and  integration  by  parts  that 

ei(t)  —  ei(0)  =  f  e^r)  dr  =  2Re  f  f  (edtU  ■  U  +  pdtV  ■  V)  clrdr 

Jo  Jo  Jn 

=  2Re  [  [  ((V  x  V)  ■  U  -  (V  x  U)  ■  V)  drdr 

Jo  Jn 

=  2Re  /  [  ((V  x  U)  ■  V  -  (V  x  U)  ■  V)  drdr  =  0. 

Jo  Jn 

Hence  we  have 


ik1/2^(-,*)iil2(n)  +  ii^1/2n,*)ii 


in) 


=  \\e1/2En 


(O) 


\\t1/2H0 


(Q)’ 


which  implies 


Ill’ll L2(Q)  +  Il''/lli2(f2)  <  H-Eo||l,2(n)  +  IIHolli2^). 
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Taking  the  first  and  second  partial  derivatives  of  (3.7)  with  respect  to  t  yields 


and 


<  dtU  +  p idfv  =  0,  V  xdtV-  edfU  =  0 

in  Cl,  t  >  0, 

:  dtU  =  0 

on  S  U  T+,  t  >  0, 

7|t=o  =  £-:(V  x  H0),  dtV\t=o  =  x  E0) 

in  Cl, 

V  x  dfU  +  iidfV  =  0,  V  x  dfV  -  ed?U  =  0 

in  Cl,  t  >  0, 

n  x  dfU  =  0 

on  su  rj,  t  >  0, 

8t2[/b  =  -M-1(Vx  (VxE0)) 

in  Cl, 

dfV\t=o  =  ~(e^)_1(V  x  (V  x  Hq)) 

in  Cl. 

Considering  the  energy  functions 

e2(t)  =  \\£1/2dtU{-,t)\\2L2{n)  +  \\p1/2dtV(-,t)\\2LHa), 

and 


es  {t)  =  \\s1/2dtU(-,t)\\2L2(n)  +  ||^t1/2i9f2V(-,t)||^2(n) 

for  the  above  two  problems,  respectively,  we  can  follow  the  same  steps  for  proving  the 
first  inequality  to  derive  the  other  two  inequalities.  □ 

4.  The  reduced  Maxwell  equations.  In  this  section,  we  derive  the  main  re¬ 
sults  of  this  work,  which  include  the  well-posedness  of  the  reduced  problem  (2.5)  and 
the  related  a  priori  estimates. 

4.1.  Well-posedness.  Let  e  =  E  U  and  h  =  H  V .  It  follows  from  (2.5) 
and  (3.7)  that  e  and  h  satisfy  the  system: 


(4.1) 


V  x  e  +  p,dth  =  0,  V  x  h  —  edte  =  J  in  Cl,  t  >  0, 
n  x  e  =  0  on  S,  i>0, 

< 

e|t=o  =  0,  h\t-Q  =  0  in  Cl, 

£?er+  =hxn+Vxn  on  T^,  t  >  0. 


Let  e  =  Jz?(e)  and  h  =  J?(h).  Taking  the  Laplace  transform  of  (4.1)  and  eliminating 
h ,  we  obtain 

!V  x  ((s/^)-1)V  x  e)  +  see  =  ~J  in  Cl, 
n  x  e  =  0  on  S, 

s-1(V  x  e)  x  n  +  3$ev+  =  V  x  n  on  T^. 

Lemma  4.1.  The  problem  (4.2)  has  a  unique  weak  solution  e  £  Hs( curl,  O), 
which  satisfies 

ll^7  x  ^lli2(n)  +  lls®lli2(n)  ~  si 1  (^11  •sJr||x.2(0)  +  \\sV  x  «||//-i/2(diVjr+) 

+  |||s|2y  x  «-||jf-i/2(div  r+)^. 
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Proof.  The  well-posedness  of  the  solution  e  £  Hg( curl,  f 1)  follows  directly  from 
Theorem  3.1.  Moreover,  we  have 

aTH(e,e)  =  (V  x  n,  er+)  —  J  ■  edr. 

R  Jn 

By  the  coercivity  of  oth  in  (3-5)  and  the  trace  theorem  in  Lemma  D.l, 

|-p-(j|V  x  e||^2(Q)  +  ||se||^2(n))  <  \\a  1  J\\l2(Q) llse||i2(0) 

+  HV'  x  p+) ||e||ij(Curi, Q) 

~  IIs  1Jr||L2(n)l|se||i,2(n)  +  || V"  x  «||jj-i/2(div  r+)||V  x  e\\L2(Q) 

+  ||s  V  X  Till  jj-l/2(d;Vj  p+)  ||sc||l,2(n)) 
which  completes  the  proof.  □ 

To  show  the  well-posedness  of  the  reduced  problem  (2.5),  we  make  the  following 
assumptions  for  the  initial  and  boundary  data: 

(4.3)  E0,H0  £  i?(curl,  Cl),  J  £  H\0,T-,  L2(L!)),  J|t=0  =  0. 

Theorem  4.2.  The  problem  has  a  unique  solution  ( E,H )  such  that 

E  £  X2(0,  T;  Li's  (curl,  ft))  n  H\0,T-,  X2(ft)), 

H  £  i2(0,  T;  H (curl,  f 2))  n  ff^O,  T;  i2(fl)). 

Moreover,  they  satisfy 

te[l?T](ll^t"^llI'2(n)  ^  X  ^'H-L2(n)  W^tH\\L2^  +  || V  x  H\\L2^) 

<  H-Eo||//(curl,n)  +  ||-H’o||jT(curl,fi)  +  II  J\\  H 1  (0,T;L2  (O))  • 

Proof.  Recall  the  decomposition  E  =  U  +  e  and  H  =  V  +  h  where  (U,  V)  satisfy 
(3.7)  and  (e,h)  satisfy  (4.1).  Since 


it  suffices  to  estimate  the  integral 

Jo  e~2sit  (llv  x  elli2(n)  +  ll^*elli2(n))  dt- 
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Taking  the  Laplace  transform  of  (4.1)  yields 

!V  x  e  +  s^h  =  0,  V  x  h  —  see  =  J  in  fl, 
n  x  e  =  0  on  S, 

3§er+  =hxn+Vxn  on  T^. 

By  Lemma  4.1, 

11^  x  e|li2(fi)  +  lls®llz,2(n)  ~  si  (lls«^lln2(n)  +  ll5^  x  nll  j/-i/2(div,r+) 

(4.5)  +  IIN  V  x  nlljf-1/2(div,r+))- 

By  (4.4), 

IV  x  ^lli2(fi)  +  lls^-lli2(n)  ~  si 1  (ll^'lli2(n)  +  lls«^llz,2(n)  +  II V  x  nllfj-i/2(div,r+) 

(4.6)  +11 M  V  x  n||Jf-i/2(diV)-]r+))- 

It  follows  from  [33,  Lemma  44.1]  that  (e,  h)  are  liolomorphic  functions  of  s  on  the 
half  plane  si  >  7  >  0,  where  7  is  any  positive  constant.  Hence  we  have  from  Lemma 
A. 2  that  the  inverse  Laplace  transform  of  e  and  h  exist  and  are  supported  in  [0, 00]. 

Denote  by  e  =  2z?_1(e)  and  h  =  2zf ~l(h ).  Since  e  =  «£?(e)  =  J^( e~Slte ),  where 
&  is  the  Fourier  transform  in  s2,  we  have  from  the  Parseval  identity  and  (4.5)  that 

/  e  2sit  (iV  x  elli2(n)  +  ll^elli2(o))  dt  =  27t  f  f||V  x  e\\2L2^  +  ||se||^2(n^  ds2 

J  O'  '  J  —00  '  ' 

/OO  pOO 

lls,^ll n2(n)^s2  +  si  /  (||sV'  x  nll jj-i/2(div, r+) 

-OO  J  —OO  '  ’  R 

+  \\\s\2V  xn\\2H^/2{d.v  T+R)yS2. 

Since  J|t=o  =  0  in  O,  V  x  n |t=0  =  dt(V  x  n) |t=0  =  0  on  we  have  J?(dtJ)  =  sJ 
in  12  and  «5f(<9t(V  x  n))  =  sV  x  n  on  It  is  easy  to  note  that 

|s|2V'  x  n  =  (2si  —  s)sV'  x  n  =  2siJf(dt,(V  x  n ))  —  ££(d2{V  x  n))  on  r^. 
Hence  we  have 

e-2sit  (||V  x  efL2{n)  +  ||ate|||2(n))  d t 

<  Sr2/_“  +  \\^id2V  x  n)HL-V2(diVir+))  ds2 

/OO 

\\J?(dtVxn)fH_1/2  +ds2. 

-OO  R 

Using  the  Parseval  identity  again  gives 

Jo  e~2sit  (||V  x  e||^2(n)  +  ||3te|||2(n))  d t 

^  «r2  /  e_2sit  {\\9tJ\\2L2(n)  +  ||0t2V  x  n||^_1/2(div  r+))  d t 
J  0 

poo 

+  (1  +  sr2)  /  e-2sit\\dtV  x  n\\2H_1/2  r+)d t, 

JO  ^  ’  R> 
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which  shows  that 


e  G  L2( 0,  T;  ffs(curl,  ft))  n  H\0,  T;  L2(fl)). 

Similarly,  we  can  show  from  (4.6)  that 

h  G  £2(0,T;if(curl,ft))n  H1  (0,T;  L2 (fl)). 

We  next  prove  the  stability.  Let  E  be  the  extension  of  E  with  respect  to  t  in 
such  that  E  =  0  outside  the  interval  [0,  t\.  By  Lemma  A.l  and  Lemma  2.2,  we  get 


Re 


—2s\t 


,rR 


&[EV+]  ■  Er+  dydt  =  Re 


lr+Jo 


e~2sit£r[Er+]  ■  Er+  dtdy 


=  —  /  Re(^[Sr+](s),£lr+(s))  ds2  >  0, 

Z7 T  J  R  R 


which  yields  after  taking  si  — >  0  that 


(4.7) 


£T[Er+]  ■  Er+  dydt  >  0. 


For  any  0  <  t  <  T,  consider  an  energy  function 

e(t)  =  \\e1/2E(;t)\\2L2{n)  +  Wfi1/2  E[(-,t)\\2L2(Qy 
It  is  easy  to  note  that 


j\\t)At=(\\e^2E{.Ml^)  +  \\^H{ 

">  ^)lli2(n) 


(4.8) 


-  (lk1/2^0|ii2(n)  +  Hm1/2^o||  i2(  a) 


On  the  other  hand,  we  have  from  (2.5)  and  (4.7)  that 


e'(t)dt  =  2Re 
=  2  Re 


iq  J a 

rt 


(edtE  -E  +  iidtH-H)  dr  elf 


o  Jn 


((V  x  H)  ■  E  -  (V  x  E)  ■  H)  drdf  -  2Re  /  /  J  -Edrdt 

Jo  Jn 

t 


=  —  2Re  /  /  &[Ev+\  ■  Er+  dydf  —  2Re 

Jo  Jr 
r* 

(4.9)  <  —  2Re 

/  o  J  n 


J  ■  E drdf 


o  Jn 


J  ■  Edrdt  <  2  max  ||-E(-,t)||i2(n)|| 


Taking  the  derivative  of  (2.5)  with  repsect  to  f,  we  know  that  {dt.E,  dtH)  satisfy 
the  same  set  of  equations  with  the  source  and  the  initial  condition  replaced  by  Dt  J, 
dtE\t=0=e-1VxE0,dtH\t=0  =  — fi  1V  x  Ha.  Hence  we  may  follow  the  same  steps 
as  above  to  obtain  (4.9)  for  (dtE,dtH),  which  completes  the  proof  after  combining 
the  above  estimates.  □ 
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4.2.  An  a  priori  estimate.  In  what  follows,  we  shall  take  a  different  route  and 
study  the  Maxwell  equations  directly  in  the  time  domain.  The  goal  is  to  derive  a 
priori  stability  estimate  for  the  electric  field  with  a  minimum  regularity  requirement 
for  the  data  and  an  explicit  dependence  on  the  time  variable. 

Eliminating  the  magnetic  field  in  (2.1)  and  using  the  transparent  boundary  con¬ 
dition  (2.24),  we  consider  the  initial-boundary  value  problem  in  a  bounded  domain: 


(4.10) 


where 


'sdfE  =  x  (/i^V  xE)-F 
E\t=o  =  Eq,  dtE\t-o  =  Ei 
n  x  E  =  0 

JVx£)  x  n  +  tf[Er+\  =  0 


in  fi,  t  >  0, 
in  Cl, 

on  S,  t  >  0, 
on  T+,  t  >  0, 


F  =  dtJ ,  Ei  =  £-1(V  x  H0  —  J0). 


The  variational  problem  is  to  find  E  G  if  5  (curl,  fi)  for  all  t  >  0  such  that 


sOfE  ■  w  dr  =  — 


(4.11) 


j  p  1  (V  x  E)  •  (V  x  w)  dr 

J  Q 

—  (rrf[Er+],w r+)  —  j  F  wdr  for  all  w  £  iTs(curl,  f2). 
R  R  J  n 


To  show  the  stability  of  its  solution,  we  follow  the  argument  in  [33]  but  with  a 
careful  study  of  the  TBC.  The  following  lemma  is  useful  for  the  subsequent  analysis. 
Lemma  4.3.  Given  f  >  0  and  any  E  £  L2(0,£;  iT_1^2(curl,  T)),  there  holds 


Re 


<^[Et+\(t)6.t\  •  Er+  (t)  dydt  >  0. 


Proof.  Let  E  be  the  extension  of  E  with  respect  to  t  in  R.  such  that  E  =  0  outside 
the  interval  [0,  £].  We  obtain  from  (A.l),  Lemma  A. 2,  and  Lemma  2.2  that 


rZ  /  rt 

—  2sit 


i  W  Jo 


<^’[Elr+](T)dT  )  •  EJr+(t)dtd7 


r+  Jo 


—2s\ t 


—2s\ t 


<^’[^r+](r)dr  )  •  Er+(t)  dtdj 


Ir+Jo 


2zf  1  o  s^o^f^r+(r)dr  •  EJr+(f)dtd7 


1 

27T  , 
1 

27T  . 


o  Er+  (s)  •  ££ (Er+ )  (s)  dyds2 

'-oo  ^r+  R  R 

(38\Ey+\,  Tlp+)ds2  >  0. 


The  proof  is  completed  by  taking  si  — ►  0  in  the  above  inequality.  □ 

Theorem  4.4.  Let  E  £  Hg( curl,  Cl)  be  the  solution  of  (4.11)  for  any  t  >  0. 
Given  Eq,Ei  £  L2(fl)  and  F  £  L1  (Q,T\  L2 (Cl))  for  any  T  >  0,  there  holds 


\E\\L°°(o,T-.L2(n))  ~  ll-®o||i2(n)  +  lli,2(n)  +  UTlIrRo ,t-l2(q.)) 


(4.12) 
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(4.13)  ||^||i,2(0,T;i2(n))  ~  T1/2||£!o||i2(n)  +  T3J2\\Ei\\L2^  +  r3/2||F||I,i(0iT;£/2(n)). 

Proof.  Let  0  <  £  <  T  and  define  an  auxiliary  function 

(4.14)  ip(r,t)  =  J  E(r,r)dr,  re!l,  0  <  t  <  £. 

It  is  clear  that 

(4.15)  V’(r  j  £)  =  0,  dtxj)(r  ,t)  =  —E(r,t). 

For  any  <j>(r,t)  €  L2(0,f-L2(il)),  we  have 

(4.16)  J  (j>(r,t)  ■  fp(r,t)dt  =  J  (^J  (j>(r,  r)dr^  •  E(r,  t)dt. 

Indeed,  using  integration  by  parts  and  (4.15),  we  have 

<j>(r,  t)  ■  tp(r,  t)dt  =  J  ^0(r,  t)  ■  J*  E{r ,  r)dr^  cl t 

=  J  Jt  E{r,r)dT-d(^J  </>(r,r)clr^ 

=  J  J*{r,T)dT.J  E(r,r)dr  +  J  <f>(r,  r)dr^  •  E(r,  t)  df 

=  Jo  (/  <^(r’T)dT)  '  ^(r^)di- 

Next,  we  take  the  test  function  w  =  ip  in  (4.11)  and  get 

f  ed2E  ■  xpdr  =  —  f  ^_1(V  x  E)  •  (V  x  xf)  dr 

(4.17) 

—  f  ^[Er+]-xjj r+ dy  —  f  F-fjjdr. 


It  follows  from  (4.15)  that 

Re  f  I  d2E  ■  xf  drdf  =  Re  f  f  ( dt(dtE  ■  fjj)  +  dtE  ■  E)  dfdr 
Jo  Jn  Jn  Jo 

=  Re/n((aE.«|)+‘i^Q<i, 

(4-18)  =l\\E(-,mh{n)-l\\Eo\\2L^n)-^  [  E1(r)  •  V>(r, 0) dr. 

J  o 

Integrating  (4.17)  from  t  =  0  to  t  =  £  and  taking  the  real  parts  yields 

l\\E(;0\\hm  -  |ll^o|||2(n)  +  V  x  E(r,t)dt\2dr 

=  eRe  f  E\{r)  ■  i/>(r,  0)  dr  —  Re  f  f  F-ff  dr 
Jo.  Jo  Jo. 

(4.19)  -  Re  t  [  <€{ Ev+ }  •  V>r+  d7df. 

Jo  J r+  R  R 
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In  what  follows,  we  estimate  the  three  terms  on  the  right  hand  side  of  (4.19)  separately. 
We  derive  from  the  Cauchy-Schwarz  inequality  that 


Re  /  E\{r)  ■  ip(r,  0)  dr  =  Re  /  Ei(r)  ■  I  /  E(r,t)dt\dr 
J  n  Wo  ' 

(4.20)  =Re  /  E\{r)  ■  E{r,  t)  drdt  <  ||-Ei|I.L2(n)  t  ||£7(-,  i)||i3(n) 
Jo  Jo  Jo 

For  0  <  t  <  £  <  T,  we  have  from  (4.16)  that 


df. 


Re 


/ o  Jn 


F1  •  if)  drdt  =  Re 
.  «  ft 

=  Re 


n  Jo  \Jo 
F(r ,  t)  •  FJ(r,  f)  clrdrdf 


JF’(r,  r)dr  )  •  E(r,  t)  dtdr 


/ 0  JO  JQ 

K  /  rt 


(4.21) 


<J^  ^||F(.,t)||ia(n)d^  ||J5(-,i)|| 

<  \\F(.,t)\\LHU)d?J  \\E{;t)\\ 


L2(i 2)uf 

i2(n)dA  • 


Using  Lemma  4.3  and  (4.16),  we  obtain 


(4.22) 


Re  /  &  (E r)  •  "i^pdydt 

Jo  Jr 

r  /•?  /  rt 


=  Re 


/rt  Jo  \Jo 


^[£ir+](r,  r)dr  I  •  F?r+ (r,  t)dMy  >  0. 


Substituting  (4.20)-(4.22)  into  (4.19),  we  have  for  any  £  G  [0,  T]  that 


dl-E(u0lli,2(fi)  +  ~ 


V  x  E(r,t)dt 


dr 


(4.23)  <  ^ll-^oll i2(0)+^y  IT(‘U)lli2(n)dt  +  e||Fii||i,2(n)^  y  ||Fi(-,t)||i2(n)dt. 

Taking  the  L  00 -norm  with  respect  to  £  on  both  sides  of  (4.23)  yields 

ll-®llit»(o,T;i2(a))  ~  ll-®o||i2(n)  +T  (||-F’||.L1(o,T;.L2(n))  +  ll-®i||i,2(n))  ll-®llz,°“(o,T;.L2(n))> 

which  gives  the  estimate  (4.12)  after  applying  the  Cauchy-Schwarz  inequality. 

Integrating  (4.23)  with  respect  to  £  from  0  to  T  and  using  the  Cauchy-Schwarz 
inequality,  we  obtain 

ll-E'lli2(0,T;i2(f2))  ~  ^ll-E'olli,2(n)+^3^2  (ITIIlRO.T;!,2^))  +  H-El||i2(n))  H-E'lll,2(0,T;L2(n))i 

which  implies  the  estimate  (4.13)  by  using  the  Cauchy-Schwarz  inequality  again.  □ 
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In  Theorem  4.4,  it  is  required  that  Eq,E\  G  and  F  €  T1(0,  T;  L2(Q)), 

which  can  be  satisfied  if  the  data  satisfy 

E0  G  L2(n),  Ho  G  H (curl,  fi),  J  G  H1  {0,T;  L2 {&,)). 

It  is  important  to  point  out  that  the  estimates  in  Theorem  4.2  were  derived  from  a 
usual  energy  method,  while  the  results  in  Theorem  4.4  were  obtained  by  using  different 
test  functions  (cf.  (4.14)). 

5.  Concluding  remarks.  In  this  paper,  we  studied  the  time-domain  Maxwell 
equations  in  a  three-dimensional  cavity.  The  open  cavity  scattering  problem  was  re¬ 
duced  to  an  initial-boundary  value  problem  in  a  bounded  domain  by  using  the  exact 
transparent  boundary  condition.  The  reduced  problem  was  shown  to  have  a  unique  so¬ 
lution  and  its  stability  was  also  presented.  The  proofs  were  based  on  the  examination 
of  the  time-harmonic  Maxwell  equations  with  complex  wavenumber  and  the  abstract 
inversion  theorem  of  the  Laplace  transform.  Moreover,  by  directly  considering  the 
variational  problem  of  the  time-domain  Maxwell  equations,  an  a  priori  estimate  was 
derived  with  an  explicit  dependence  on  time  for  the  electric  field.  Computationally, 
the  variational  approach  leads  naturally  to  a  class  of  finite  element  methods.  As  a 
time-dependent  problem,  a  fast  and  accurate  marching  technique  shall  be  developed 
to  deal  with  the  temporal  convolution  in  the  transparent  boundary  condition.  We  will 
report  the  work  on  its  numerical  analysis  and  computation  in  a  forthcoming  paper. 

Appendix  A.  Laplace  transform. 

For  any  s  =  Si  +  is2  with  si,  S2  G  R,  Si  >  0,  define  by  u(s)  the  Laplace  transform 
of  the  vector  field  u(t),  i.e., 

/>00 

it(s)  =  Jf(u)(s)  =  /  e~stu(t)dt. 

Jo 

It  can  be  verified  from  the  integration  by  parts  that 

(A.l)  [  u(t)  dr  =  T£~x 

Jo 

where  Jz?-1  is  the  inverse  Laplace  transform, 
inverse  Laplace  transform  that 

(A. 2)  u(t)  =  &~1(eSltJ£{u)(s\  +  is2)) , 

where  ^_1  denotes  the  inverse  Fourier  transform  with  respect  to  S2- 

Recall  the  Planclrerel  or  Parseval  identity  for  the  Laplace  tranform  (cf.,  [13, 
(2.46)]). 

Lemma  A.l.  Ifu=J£(u)  andv  =  J£(v),  then 

pOO 

ii(s)  ■  u(s)ds2  =  /  e~2situ(t)  ■  v(t)dt, 

Jo 

for  all  si  >  A,  where  A  is  the  absissa  of  convergence  for  the  Laplace  transform  of  u 
and  v. 

The  following  theorem  (cf.,  [33,  Thm  43.1])  is  an  analog  of  the  Paley-Wiener- 
Schwarz  theorem  for  the  Fourier  transform  of  the  distributions  with  compact  support 
in  the  case  of  Laplace  transform. 

Lemma  A. 2.  Let  w(s)  denote  a  holomorphic  function  in  the  half  plane  si  >  <7o, 
valued  in  the  Banach  space  E.  The  following  statements  are  equivalent: 


(s  ^(s))  , 

One  verifies  from  the  formula  of  the 
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1.  there  is  a  distribution  w  £  V'+{E)  whose  Laplace  transform  is  equal  to  w(s); 

2.  there  is  a  <j\  with  <Jq  <  a\  <  oo  and  an  interger  m  >  0  such  that  for  all 
complex  numbers  s  with  si  >  <j\,  it  holds  ||'Uj(s)||e  <  (1  +  |s|)m, 

where  V'+{E)  is  the  space  of  distributions  on  the  real  line  which  vanish  identically  in 
the  open  negative  half  line. 

Appendix  B.  Spherical  harmonics  on  hemisphere. 

The  spherical  coordinates  (p,  6 ,  y>)  are  related  to  the  Cartesian  coordinates  r  = 
{xi,X2,xs)  by  X\  =  p  sin  f?  cosy),  X2  =  psm.9 sirup,  X3  =  pcosO,  with  the  local  or¬ 
thonormal  basis  {ep,  eg,  ev}~. 

ep  =  (sin 0 cosy;,  sin 8 sin <p,  cos 6), 
eg  =  (cos  8  cosy?,  cos#  sin  y>,  —  sin#), 
ev  =  (- siny>,  cosy),  0), 

where  8  and  p  are  the  Euler  angles.  Let  T  =  {r  :  p  =  1},  T+  =  {r  :  p  =  1,  X3  > 
0},  T“  =  {r  :  p  =  1,  X3  <  0}  be  the  unit  sphere,  upper  unit  hemisphere,  and 
lower  unit  hemisphere,  respectively.  Denote  by  T p>  =  {r  :  p  =  R},  =  {r  :  p  = 

R,  X3  >0},  T~  =  {r  :  p  =  R,  x 3  <  0}  the  whole  sphere,  upper  hemisphere,  and  lower 
hemisphere  with  radius  R,  respectively. 

Let  {FJ71(0,  p),  \m\  <  n,  n  =  0, 1, 2, . . .  }  be  an  orthonormal  sequence  of  spherical 
harmonics  of  order  n  on  the  unit  sphere  T  that  satisfies 

(B.l)  ArY™  +  n(n+l)Y™  =  0, 

where 

.  1  d  (  .  n  d  \  1  d2 

r  “  sin  8  88  \Sm  °  88  )  +  sin2  8  dp 2 

is  the  Laplace-Beltrami  operator  on  T.  Explicitly,  the  spherical  harmonics  of  order  n 
is  written  as 


Y™{8,p)  = 


/2n  +  1  (n—  |m|)! 


47T 


|m|)! 


Plm\cos9)ei‘ 


where  the  associated  Legendre  functions  are 

Amp  (f\ 

Pn(t)  :=  V(1  -t2)m  m  =  0,1, . . .  ,n. 

Here  Pn  is  the  Legendre  polynomial  of  degree  n,  which  is  an  even  function  if  n  is  even 
and  an  odd  function  if  n  is  odd. 

Define  a  sequence  of  rescaled  spherical  harmonics  of  order  n : 

(B.2)  X™{9,p)=V^Y™{9,p), 

which  forms  a  complete  orthonormal  system  in  L2(T^)  for  \m\  <  n,  m  +  n  =  odd,  n  £ 
N  (cf.,  [29,  Lemma  3.1]). 


20 


P.  Li,  L.  Wang,  and  A.  Wood 


Denote  by  L2(T^)  the  complex  square  integrable  functions  on  the  hemisphere 
For  convenience,  we  take  the  following  notation  for  double  summations: 


oo  n 


E  <  ■  =  E  E  < 

\m\<n  n—1  m=—n 

odd  oo  n 

E  <■ ;  =  E  E  ■ 

|m|<n  n—  1  rn=—n 

m+n— odd 


E  :  =  E  E 

|m|<n 


n—1  rn——n 
m+n— even 

To  describe  vector  wave  functions  on  the  hemisphere,  we  introduce  some  boundary 
differential  operators.  For  a  smooth  scalar  function  w  defined  on  1+  let 

dw  1  dw 

vrw  =  -xw  eg  + 

off  sin  6  dip 

be  the  tangential  gradient  on  The  surface  vector  curl  is  defined  by 

curlrw  =  Vr w  x  ep. 

Denote  by  divp  and  curlp  the  surface  divergence  and  the  surface  scalar  curl,  respec¬ 
tively.  For  a  smooth  vector  function  w  tangential  to  1+  it  can  be  represented  by  its 
coordinates  in  the  local  orthonormal  basis: 


w  =  we  eg  +  wv  e<p, 


where 


wg  =  w  ■  eg  and  wv  =  w  ■  ev. 

The  surface  divergence  and  the  surface  scalar  curl  can  be  defined  as 

1 


divpio  = 
curlp  in 


- :  -  0 
1 

sind 


9  ,  ■  as  ,  9wv 

9  dwe 


It  is  known  (cf.,  [32])  that  these  boundary  differential  operators  satisfy 

(B.3)  Ap  =  divpVr  = —curlp  curlp  and  curlp  Vr  =  divpcurlp  =  0. 

It  is  also  know  (cf.,  [14,  Thm  6.23])  that  an  orthonormal  basis  for  X2(r+  =  {in  £ 
L2(Tfi)  :  ep  ■  w  =  0},  the  tangential  fields  on  r^,  consists  of  functions  of  the  form 

1 


Ry/n(n  +  1) 


vi+r++ 


and 


V™  (ff,  ip)  =  ep  x  U™  (9,  ip)  =  — 


R^/n(n  +  1) 


curlvY. * 
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for  |m|  <  n,  n  €  N.  It  follows  from  (B.l)  and  (B.3)  that 


and 


divrtC 


\Jn(n  +  1) 
R 


ym 

1  n  i 


CUrlrV™  = 


sjn(n  +  1) 


ym 

1  n  i 


curlr  U™  =  divpV”  =  0. 

Define  two  sequences  of  tangential  fields 

(B.4)  X™(9,  <p)  =  1  VrJTO  ip)  =  V2  U™(9,  ip) 

V  nyn  +  1) 

and 


Y™(0,  p)=epx  X™(0,  <p)  =  V2  V™(0,  p). 


Using  the  definition  of  the  tangential  gradient,  and  noticing  that  eg  x  ev  =  epi  ev  x 
ep  =  eg ,  ep  x  eg  =  ev,  we  get 

7T 

eg  x  X™(—,p)  =  0  for  |m|  <  n,  m  +  n  =  odd,  n  €  N, 


and 


7 r 

eg  x  Y™(—,tp)  =  0  for  \m\  <n,m  +  n  =  even,  n  €  N. 

Define  a  subspace  of  complex  square  integrable  tangential  fields  functions  on  the 
hemisphere  rj  : 

Lt(TR)  =  {w  e  i2(r+) :  ep  ■  w  =  0}. 


It  is  shown  (cf.,  [29,  Lemma  3.2])  that  the  vector  spherical  harmonics  { X ™  :  m  +  n  = 
odd}  and  {Y™  :  m  +  n  =  even}  for  |m|  <  n,  n  €  N  form  a  complete  orthonormal 
system  in  X^(r^). 

Appendix  C.  Hankel  functions. 

For  v  £  R,  the  two  Hankel  functions  H^\z)  and  Hl2\z)  where  z  €  C,  are  two 
fundamental  solutions  of  the  Bessel  equation  for  functions  of  order  v\. 


2d  2u  dti  2  2n 

z  -TT  +  Z—  +  {z2  -  v2)u  =  0. 
dzz  dz 


Recall  the  Bessel  funtions  of  imaginary  argument  Kv(z),  also  called  the  modified 
Bessel  functions,  which  is  the  solution  of  the  differential  equation 


dii 


Z~ -pw  +  Z- - ( Z 2  +  V2)u  =  0. 


dz2  dz 

It  is  connected  with  H„\z)  through  the  relation 


(C.l) 


Kv{z)  =  -Tri  e^H^iiz). 
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It  is  known  (cf.,  [38,  p.  511])  that  Kv(z)  has  no  zeros  if  |argz|  <  *7r,  which  implies 
from  (C.l)  that  H^1\z)  has  no  zeros  when  Imz  <  0. 

The  spherical  Hankel  function  hn^  (z)  can  also  be  defined  by  the  Hankel  function 
of  half  integer  order: 

(c.2)  '■")<2)  = 

Comibining  (C.l)  and  (C.2)  yields 

h£\z)  =  -y^i 

which  implies  that  hn\z)  has  no  zeros  when  Imz  >  0. 

The  following  two  lemmas  on  the  spherical  Hankel  functions  for  the  complex 
number  are  proved  in  [11]. 

Lemma  C.l.  Let  R  >  0,  n  £  Z,  s  =  si  +  is2  with  s\  >  0.  R  holds  that 

Re{  1  +  (i sR) }  <  0, 

where  r„\z )  =  z(hn))'(z)/hn\z). 

Lemma  C.2.  Let  R  >  0,  n  £  Z,  s  =  si  +  is2  with  si  >  0.  Then  we  have 
Im{s2(l  +  r^(is.R))}  <  0. 


Combining  the  result  in  Lemma  C.l  for  small  value  of  n  and  the  proof  in  [14] 
(cf.  [25,  Lemma  3.1])  for  large  value  of  n,  we  may  obtain  the  following  estimate. 

Lemma  C.3.  Let  R  >  0,  n  €  Z,  s  =  s±  +  is2  with  si  >  0.  There  exist  positive 
constants  C\  and  Ci  such  that 

Cm  <  |1  +  (isi?)  |  <  Cm- 


Appendix  D.  Functional  spaces. 

Denote  by  1LS(T^)  the  Sobolev  space,  the  completion  of  C“(T^)  in  the  norm 
II  ‘  characterized  by 

odd 

IHIffs(r+)  =  E  (1  +  n(n  +  l))s  \w™\2, 

|m|<n 

where 

odd 

w(0,<p)=  E  <**n  (*>¥>)■ 

|m|  <n 

Introduce  three  tangential  trace  spaces: 

Hst(T r)  =  £  (HS{T+))3,  ep  ■  w  =  0,  eg  x  w(tt/2 ,ip)  =  0}, 
ff“1/2(curl,  T+)  =  {wG  H~1/2(T+),  curlrtc  £  LT1/2(T+)}, 

H~1/2( div,  T+)  =  {«,  £  h;1/2(t+),  divrw  £  H -1/2(r+)}. 
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For  any  tangential  field  w  G  it  can  be  represented  in  the  series  expansion 


odd  even 

t»=  E  <nX^(9,<p)+  Y  w?nYZ(0,<p). 

|m|  <n  |m|<n 

Using  the  series  coefficients,  the  norm  of  the  space  ff®( r^)  can  be  characterized  by 

odd  even 

IHlLt*(r+)  =  (!  +  n{n+  i))s  \w^J2  +  Y  (1  +  n{n  +  l))s  |w^|2; 

|m|<n  |m|  <n 

the  norm  of  the  space  i?_1^2( curl,  T^)  can  be  characterized  by 

odd  1  even 

ll“ll»-/.(a,,.,ri)  =  E  7TT  7  J.,.K‘nl2  +  E  ✓l  +  n(n  +  l)K”|2; 

|m|<n  V  '  /  |m| <n 

and  the  norm  of  the  space  i/-1/2(div,  T^)  can  be  characterized  by 

odd  even 


lif-1/2(div,r+) 


=  Y  V1  +  n{n+l)\w^n\2  +  Y 


1 


|m|<n 


|m|<n 


y/l  +  n(?l  +  1) 


Kl2- 


Define  a  dual  pairing  by 


where 


(u,v) 


+  Y  U™nV  Tn 

|m|  <n 


odd  even  odd  even 

«  =  Y  UTn*n+  Y  <nYn  and  «  =  Y  <n  Xn  +  Y  V&YZ ' 

|  m|  <n  |m|<n  |m|<n  |m|  <n 


Introduce  three  functional  spaces 

H  (curl,  Cl)  =  {u  £  r2(U),  Vx£e  £2(U)}, 
i?s(curl,  f2)  =  {it  €  iT(curl,  fi),  n  x  u  =  0  on  S'}, 
ffo(curl,  Cl)  =  {u  £  H( curl,  fi),  n  x  «  =  0  on  SUT^}, 


which  are  Sobolev  spaces  with  the  norm: 

1 1 u  1 1 H  (curl,  Q)  =  ^ll'ulll,2(n) 

Given  a  vector  field  u  on  r}j,  denote  by 


itr+  =  —ep  x  ( ep  x  u) 

the  tangential  component  of  u  on  which  satisfies  the  following  trace  estimate 
(cf.  [29,  Lemma  3.3]). 

Lemma  D.l.  For  any  u  £  iTs(curl,  Cl),  there  holds 


llMr+ll j/-1/2(curi,r+)  ~  llMll«(curi,n)- 
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